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Abstract. Given a function ip in C 2 (M. d ), the affine (wavelet) system gener- 
ated by rjj, associated to an invertible matrix a and a lattice T, is the collection 
of functions {| det a\i l 2 i\>{a? x — 7) : j £ Z, 7 G T}. In this article we prove that 
the set of functions generating affine systems that arc a Riesz basis of C 2 (M. d ) 
is dense in £ 2 (R d ). 

Wc also prove that a stronger result is true for affine systems that are a 
frame of £ 2 (R d ). In this case we show that the generators associated to a fixed 
but arbitrary dilation are a dense set. 

Furthermore, we analyze the orthogonal case in which we prove that the 
set of generators of orthogonal (not necessarily complete) affine systems, that 
are compactly supported in frequency, are dense in the unit sphere of C 2 (M. d ) 
with the induced metric. As a byproduct we introduce the p-Grammian of 
a function and prove a convergence result of this Grammian as a function of 
the lattice. This result gives insight in the problem of oversampling of affine 
systems. 



1. Introduction 

Let ip be an £ 2 (R d )-function, a an invertible d x d matrix and T a lattice in K d . 
The wavelet system or affine system (jp,a,T), generated by ip and associated to a 
and r, is the collection of functions, 

(i>,a,T) = {D aj T^:jeZ, 7 GT}. (1) 
Here T y and D a denote the the unitary operators in £ 2 (R d ) defined by (T y iji){x) = 
*l>(x -y),y£ R d and (D a ip)(x) = | det a\ 1/2 ip(ax), a G GL d (R). 

Affine systems have been studied in depth during the last 25 years mainly because 
of their importance in applications. In addition they proved to be very useful in a 
variety of theoretical problems. On the other hand they were studied in the context 
of Hilbert spaces where the translation and dilation operators were replaced by 
a general group of unitary operators [GLT93I IUL98j . Generalizations of wavelet 
systems in £ 2 (K d ) with translations not necessarily on a lattice and using different 
dilations were also considered |ACMn4| . |CLL+n4) . 

One of the relevant questions about these systems is whether the collection 
(ijj, a, T) form an orthonormal basis, a Riesz basis or a frame of the space £ 2 (R d ). 

The construction of affine systems with specific prescribed properties is a dif- 
ficult problem. Usually it is accomplished imposing conditions on the generators. 
This has been one of the core problem from the beginning of wavelet research 
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|Dau88l |Mey88[ |Mey92| IBF94I ICHM04j (and references therein). Still there are 
many unanswered questions and open problems in the study of affine systems. 

One way to get a better understanding of these systems is considering the set 
of generators of affine systems having some particular structure, and trying to 
answer global questions about these sets. For example, wavelet systems forming a 
tight frame of the space have been completely characterized by a set of equations 
imposing conditions on the Fourier transform of the generators [CCMW021 ICS00I 
IBow03| . Another example is the problem of connectivity of the set of generators 
which has received considerable attention lately, |Con 98 Spe99|. 

In this article, motivated by a question posed by David Larson |Lar05| . we study 
the problem of density of the set of generators of affine systems. We prove (Thco- 
rem l4.1[l that the set of generators forming Riesz bases is dense in C 2 (M. d ) when we 
allow diagonal matrices and arbitrary lattices. This theorem gives an even stronger 
result, since the generators in the dense set have Fourier transform supported on 
wavelet-sets and as a consequence have orthogonal dilations. The proof is obtained 
as a combination of the techniques used in the proof of density for wavelet frames 
(section |2J) and the theory of wavelets sets that has been developed recently by 
different groups of rese arches ^U)J\ IBMM99I lBSn2l lBS04l IDLS97I IDLS98I ITLP98I 
ISW98llWlm02llZak9fi| . 

The question about density of generators of wavelet frames, is also answered 
positively. We obtain a very general result, that confirms the flexibility in their 
construction. If a dilation matrix a is chosen arbitrarily, the set of frame generators 
associated to the dilation a is dense in £ 2 (M. d ) (Theorem The main tool here 
is the use of the general scheme in frame construction that appears in |ACM04j . 

Finally in section [S] we study the case of orthogonal affine systems. Since gener- 
ators of orthonormal systems have norm one the question here is whether they are 
dense in the unit sphere of C 2 (M. d ). It is easy to see that this is not the case when 
cither the dilation or the lattice is fixed (see section EJ). 

We first show that for a given function of norm one, we can always find a lattice 
r and a function with orthonormal translates in that lattice, that is close to the 
original function. We then prove that a dilation can be chosen in such a way that 
the associated affine system is orthonormal. This proves the density of generators of 
orthonormal wavelet systems (not necessarily complete in C 2 (M. d )) fTheorcm l5.2J) . 

The strategy here is the study of the behavior of the Grammian of the generator, 
as a function of the lattice. We obtain an interesting result on the convergence of 
the p-Grammian. We prove that the p-Grammian of a function converges point- 
wise and in £ p ([0, l] d ) to the constant function 1 when the lattice expands. This 
result is a statement about oversampling. 

However the question whether complete orthonormal affine systems are dense in 
the sphere remains unanswered. 



We will denote by n the Lebesgue measure in M. d and by GLd(M.) the usual group 
of invertible matrices in W lxd . By a lattice T we mean T = cE d , where c S GLd(M.). 
The Fourier transform of a function / <G C 2 (M. d ) is 



2. Notation 




dx. 
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For a measurable Q C K d , Kq will denote the functions in £ 2 (R d ) with support in 

K n = {/ G £ 2 (K d ) : supp(/) C fi}. (2) 

We will say that a matrix a is expansive if a £ GLd(M) and |A| > 1 for all eigenvalues 
A of a. 

Since affine systems depend on the matrix and the lattice involved we consider 
the following sets of generators: 

Wo{a,T) := {ip £ C 2 (R d ) : (i(j,a,T) is an orthonormal basis (onb) of C 2 (R d )} 
and 

Wo(a):= (J W (a,cZ d ) 

ceGL d (J&) 

Wo(T):= (J W (a,r) 

aeGL d (K) 

W G := (J W (a,cZ d ). 

a,ceGL d (R) 

Similarly, we use the notations Wr or ~Wp for the corresponding generators of 
wavelet Riesz basis and frames respectively. 

3. Density of the set of frame generators 

In this section we will prove that for a fixed expansive matrix a, the set VVf(i) 
is dense in C 2 (R d ). 

In what follows diam(5') will denote the diameter of a set S C M. d and B(Q,t) 
the ball in K d centered at and with radius t. 

Let X be a discrete set in M. d . The gap p of X is defined as: 

p = p( x ) = su p in t \ x ~ tI- 

The set X is said to be separated if inf 7 ^ 7 / I7 — 7' > 0. 

For a separated set X Beurling |Beu66| proved the following result: 

Theorem (Beurling). Let X C M d 6e separated, and D, = 5(0, r). //r/j < 1/4, 
t/ien {e _2lr?7 ' i 'xa(w) : 7 £ X} is a frame for /Co- 

For a very clear exposition of some of the Beurling density results see |BW99j . 

We need now the following result which is a particular case of the more general 
theorems in |AUM04j . We include its proof here for completeness of this presenta- 
tion. 

Theorem 3.1. Let a be an expansive matrix and h a bounded compactly supported 
function such that \h(uj)\ > c > 0, for almost every cu G hi := supp{h), and ^ 
supp(h). Lf X C K d is a separated set such that the gap p(X) < ^ where R is such 
that U C B(0, R) then 

{D a3 T^:j eZ,jeX} 
is a frame of C 2 (M. d ) where ifj = h. 
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Proof. By Bcurling's theorem we know that 

is a frame of Kb(q,r) with some frame bounds < m < AI < +00. In particular 

{e- 2 ^Xu(u) : 7 G X} 

is a frame of ify. Using that D a is a unitary operator and calling b = (a -1 )*, we 
conclude that for each j G Z 

{| dct(6)p/ 2 e- 2 ^^ Xw (6^) : 7 G X} 

is a frame of K^-iy with the same bounds m, M. 

Further, by Proposition 5.9 in |ACM04j there exist < p < P < +00 such that 

p < E IM^MI 2 < P, a - e - w G M d . 

Now, since supp(/i) = U, if we set hj{w) = h(b>oS) then supp(/ij) = b~ 3 U. 
Given now / £ £ 2 (R rf ) and calling fj = /ij/ we will see that 

P \\f\\ 2 < Eii/^^imi 2 , (3) 

3 

for 



p||/ 



j P \f\ 2 <Jy: i^fi/1 2 = e/ i^/i 2 = e ii/.-ii 2 



j j j 

and the upper inequality can be obtained similarly. Using that fj € K^jjj we have 

m||/#< ^\<f j (u),\det(bW/ 2 e- 2 « ibj ^ Xb -iu(")>\ 2 < M||/,-|| 2 . (4) 

7 

In addition, since supp(/j) = b~ J lA, we have that 

< /,-(«), I det^Jl^e-^^Xi-^H >=< /(o;),|det(&)[ 3 VV W ^/i>) > . 

(5) 

Using (J3J) in Q and summing in j by J2J we get 

pm||/|| 2 < ^^|</H,|det(6)p/ 2 e - 2 - b3 ^H>| 2 < PA/II/II 2 . 

3 z 9 

The claim is now a consequence of Plancherel's Theorem. □ 
Now we are ready to state and proof the main result of this section. 

Theorem 3.2. The set Wjr(a) is dense in £ 2 (M d ), for every dxd expansive matrix 
a. 

Proof. Assume that a function / S C 2 (R d ) and a positive number e are given. 
We want to approximate / by a wavelet frame function tp. We will do this by 
constructing the Fourier transform of ip. 

• Select g £ £ 2 (R d ) such that g is a continuous function and ||/ — g\\% < e/2. 
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• Choose < r < R < oo such that 

aB(0, r) C B(0, R) and / \g(u)\ 2 duj < — , 

jR d \U 8 

where U = {ui : r < \u>\ < R}. 

e 2 

• Choose A > 0, such that 4A 2 /i(W) < — , and let 

8 

E\ := {uj eR d : \g(to)\ > A}. (6) 

• Define a function h in the following way: 

!g(u) lu e U n Ex 
A u £ U \ E x ■ (7) 
else 

• Choose X to be any separated set with gap p{X) < j^. (Note that X can 
be chosen to be a lattice.) 

We will now see that we are indeed under the hypothesis of Theorem 13. II Since 

Q := aB(0, r) \ 5(0, r) C B(0, R) \ B(0, r) = U 

we note that by Lemma 5.11 of |ACM04j {a j Q} is a covering of R d \ {0} and 
therefore {aPU} is also a covering of K d \ {0}. 

Finally, h and X were defined to match the hypothesis of Theorem l3.ll Therefore 
we can apply that Theorem to conclude that (ip, a, X) is a frame for £ 2 (M. d ), with 

ip = h. 

Furthermore, 

U-h\\l= I \g-h\ 2 + { \ g -h\ 2 

JU JR d \U 




<4X 2 f x(U\E x ) + j< E — 

where the last inequality comes from the choice of A. 
Finally 

||/-^||2<||/-5l|a+||§-h||2<e. 

□ 

4. The Riesz basis case 

In this section we will prove that the set of generators for Riesz basis is dense 
when we allow diagonal matrices and arbitrary lattices. Since the requirement of 
building a Riesz basis is much stronger than to form a frame, the construction of 
the function requires more subtle techniques than for the previous case. 

The idea is again to approximate the given function on the Fourier side with 
an appropriate function. We will use the construction for the frame wavelets and 
adapt it to the Riesz basis case. The main difficulty here lies in the fact that the 
elements of the affinc system need to be independent, which forces the construction 
to be more involved. 
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It is worth to remark here, that the construction in the previous section was very 
general and stable in the sense that we had freedom to move things a little bit - 
and still obtain a satisfactory result. In this section, the construction is tight and 
very adapted to the function we want to approximate. 

The method we use, is based on a standard wavelet-set construction. Wavelet 
sets have been studied and developed by many groups (for references see the In- 
troduction). We adapt the construction of |Zak96| and |BS02) to obtain functions 
supported on a wavelet set, but which are not constant. In this way we obtain a 
Riesz basis, instead of a wavelet basis. 

In what follows we will say that two sets A and B in K rf are almost disjoint if 
(i(AnB) = 0. 

Lemma 4.1. Let ft C K d be a set of finite measure. If {\k}k£Z C R d satisfies 
that {e 2 ^ lXkU x^( UJ ) : fc G Z} is a Riesz basis for /Co with bounds A and B and 
h G C 2 (R d ) satisfies that < p < \h(cu)\ < P < +oo then 

{h{cj)e 2 " XkUJ : k G Z} 

is a Riesz basis for /Co, with Riesz bounds pAfx(fl) and PB(i(£l). 

If a G GL c i(W.) and f2 satisfies that Ujgza J 'f2 = M. d up to a set of zero measure, 
with the union being almost disjoint, and {gk '■ k G Z} is a Riesz basis for /Co, then 

{Dig k :k,j£Z} 
is a Riesz basis for C 2 (M. d ) with the same bounds. 

Proof. The first assertion is immediate, and the second one follows from the fact 
that the dilation is a unitary operator in £ 2 (K d ). □ 

We are now ready to state the main theorem of this section. 

Theorem 4.1. The set W R is dense in £ 2 (R d ). Precisely, if f G £ 2 (tt d ), and 
£ > 0, there exists a function ip G C 2 (M. d ), an expansive matrix a G R dxd , and a 
lattice T, such that 

• 11/ - i>h < £ and 

• (*,a,r) is a Riesz basis for £ 2 (R d ). 

Before proceeding with the proof let us give the following definition. 

Definition 4.1. Let T be an arbitrary lattice. A set E C M. d is T-congruent to a 
set E C R d , if there exist partitions {E s : s G T} of E and {E s : s G T} of E in 
measurable sets such that for every s G T, E s = E s + s. 

In the proof of Theorem l4.1l we will use the £°°-norm of M. d . In this way we will 
obtain an orthogonal basis of exponentials supported on the cube B 00 {Q,R/2) for 
some appropiatc R. If we multiply the elements in this basis by a function that 
is bounded above and bounded away from zero, by the previous Lemma, we will 
have a Riesz basis on /C^ 00 (o,_r/2)- Therefore, looking at the previous construction 
in the frame case, we would now need a set U on which \g\ is bounded away from 
zero, that is T-congruent to B^O, i?/2), and furthermore, that tiles the plane by 
dilations by a matrix a. This forces us to be more careful in the choice of r, and 
will also limit our choices of L and a. 
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4.1. Proof of Theorem 14.11 

Proof. Let as before g £ C 2 (M. d ) be such that g is continuous and ||/ — 5II2 < f ■ 
We will approximate g with an appropriate function. 

• Let R > be such that J Rti y B ^ R , 2 ^ \g(uj)\ 2 duj < j^. 

• We now select r > small enough such that: 

r<f (8) 



L 



\g(co)\ 2 <ko < - (9) 

3^(0,7-/2) 

e 2 

and mr d < — where m := max{|<?(tj)| 2 : lo G i?oo(0, i?/2)}. (10) 
16 

Let now 

R 

T = R1 d and a = — I dxd . 

r 

We define 

f:R d — > R d 

x 1 ► x ~\~ R£j if x G j — quadrant 

where ^ is the vertex of the cube [—1,1]" that lies in the j th -quadrant. 
Let us call 

A) = Boo(0,r/2) and define A t := (a" 1 o f Y(A ), i = 1, 2, . . . (11) 

It will be convenient to use the notation A 3 for the intersection of the set Aq with 
the j th -quadrant. With this notation, note that 

l<j<2 d 

= |J r Jlo -i C , + ... Al -^(o- i (4)). (12) 

l<j<2 d 

Here T y denotes the usual translation by y in £ 2 (IR d ). Therefore we have that 

OO OO j ■ 

/ ip \ ai f u 

EM^) = ^o)Ey = rd jn< rd > ( 13 ) 

i=l i=l 

where the last inequality holds by JSJ. 
We define the set 

and for A = 8 ^ d / 2 , the function /i by 

{<?(w) x € U C] E\ 
A x£U\E x 
else 

where as before 

£ A :={^eR d :|gH|>A}. 
Our claim is that ip with ip = h satisfies both conditions of the Theorem. 
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For the first, let us compute 

11/ - </>||2 = 11/ - h\\ 2 < 11/ - g\\2 + - h\\ 2 < \ + ||S - h\\ 2 . 
To compute \\g — h\\2 we note that \\g — h\\^ can be split into 3 integrals 

f \g-h\ 2 + [ \g~h\ 2 +[ \g-h\ 2 . (14) 

JUnE x JU\E X JR d \U 

By the definition of h, the first term in (f 1 4|) vanishes. 

For the second, note that onW\ E\, \g(u)\ < A and therefore 

I \g-h\ 2 <4\ 2 n(U\E x )<^-. (15) 
Ju\e x 16 

For the last term, since h(u>) = if oj U we have 

l5| 2 </ l5| 2 +/ l3| 2 - (16) 

JB^(0,R/2)\U JM. d \B aa (0,R/2) 

The right hand side of Q16f) can be written, 

l3| 2 - (17) 







f 


/ l.9| 2 + 


[ \9\ 2 + 






Ju i>1 A i 


■'R< 1 \B oo (0,fl/2 



The first and last term in this sum arc each smaller than e 2 /16 by our choice of r 
and R. For the middle one, we use the computation about the measure of U,>ij4j 
done in 1)1 3|) and the choice of r in (|10)l to obtain 

r £ 2 

/ \g\ 2 duj < mM(U i >iA i ) < mr d < — . (18) 
Ju.mti, 16 

Putting all this together, we obtain 

£ 2 

— ^lll < -j and therefore ||/ - -0 II 2 < £■ (19) 

This proves that the function h can be chosen as close to g as we wish. It remains 
to show, that (ip,a,T) is a a Ricsz basis. For this, we observe that: 

• By construction, U tiles M. d \ by dilations by a. 

• Furthermore U tiles R d by translations on T. For this, we first note that if 
x G A n then 

rR / /r\"\ .. rR / / r \ n ( r + R\\ 

<i - U ))< Moo < —— i - - hnr ■ (20) 



R-r\ \RJ /-" " R-r\ \RJ \ 2R 

This fact allows us to conclude that: 

(1) At C {B oo (0,R/2)\A ),i>l, 

(2) AinAj=<Di^ j 
which allows us to rewrite U 

(oo \ / oo > 

aAo\(\jAi) U MJfAi 
i=0 / \i=0 / 

This shows that U is T-congruent to Bqo(0, -R/2). 
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On the other hand, since A < < m on U by Lemma PTD 

is a Riesz basis for JCu- 

Thus we found a Riesz basis for JCu-, and IA is a wavelet set for the dilation a and 
translation T. Therefore, the system (ip,a, T) is a Riesz basis of £ 2 (R d ). 

□ 

Remark. In the above proof, we constructed a function that is supported on a 
wavelet set. In this way we in fact prove that a proper subset of Wr is dense in 
£ 2 (R d ), since the functions in the dense set are supported on a wavelet set. Further, 
they generate quasi-orthogonal affine systems, since they have orthogonal dilations 
(i.e. < D a3 T^,D aJ ,T Y ^>= if j ? f). 

5. Density and Orthogonal Wavelets 

In this section we will consider the problem of the density for orthonormal 
wavelets. Since orthonormal wavelet functions have norm one, the natural question 
in this context is wether they are dense on the unit sphere of £ 2 (M. d ) with the 
induced metric (i.e. S d := {/ G £ 2 (R d ) : |]/]| 2 = 1}). 

An immediate argument, that we will see later, shows that if we fix either the 
dilation matrix a or the lattice T, then the sets Wo (a) and Wo CO are not dense 
in Sd- If we allow both, the matrix a and the lattice T to be arbitrary, it is an open 
problem if the set Wo is dense in Sd- 

In this section we will prove that the set of functions 

Wo := {ip ■ 3 r and a G GLd(R), such that (ip, a, T) is an orthonormal system} 

is dense in Sd ■ Note that in Wo we removed the completeness requirement of the 
system. 

In particular we prove an interesting property of the Grammian of an arbitrary 
function in £ p ([0, l] d ) which is of independent interest. This property is then used 
to derive some consequences which in particular imply the above mentioned result. 

Throughout this section, b will be a matrix in GLd(R). Our first Lemma is a 
known result (see for example |Mal89j ) . We state it here in the form we need it. 

Lemma 5.1. Let f G £ 2 (R d ). The following statements are equivalent 

(1) The system {T(( h »)-i fc )/ : k E Z d } is orthonormal in £ 2 (M. d ). 

(2) The system {/(w)^ 2 "^ 6 *^^ : k G Z d } is orthonormal m £ 2 (R d ). 

(3) Zsez«\\ dctb \ 1/2 f( b (" + s )\ 2 = 1 a - e - m IM d - 

(4) Es^l/^ + MPHdetbr 1 a.e. m b([Q, l] d ). 

Proof. We will show Q ■ The rest is trivial. 

The system {r ((b .)-i fe) / : k G Z d } is orthonormal in £ 2 (M. d ), if and only if 

8o.k= f f(x-(b*)- 1 k)7(x)dx= [ e- 2 — ^ lk f\co)J{u)duj 

= [ (| det 6| V|/>(w + s))| 2 ) e~ 2 ^- k dw. 
J[o,i] d 

Hence, it will be orthonormal if and only if © is satisfied. □ 
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Note that the sum in © can be rewritten as X)sgz d \TsDbf(uj)\ 2 , which motivates 
the following definition. 

Definition 5.1. f f g £ £ 2 (R d ) and c S GLd(M.), the Grammian of g with respect 
to c is the function 

ffcM = ( ^ |T s ^ c5 H| 2 ) 1/2 = [detc| 1 /2( ^ |. ( c(w + s ))| 2 ) 1/2 

Note that g c is a Z d -periodic function. 

Tn what follows the set E(g,b) := {lo 6 M. d : <7&(w) > 0} will be relevant. It is 
immediate to verify that, if g G C 2 (M. d ), then 

Il5b||£ 2 ([0,l] d ) = ll.9ll£ 2 (R d )- ( 21 ) 

For g e £ 2 (R d ), we define u g by 



g(") 



if w e b) 



Ug (u) := I 9 b {b-^) ----- w >-v (22) 
I otherwise. 

We have the following Lemma, whose proof is immediate. 

Lemma 5.2. Let g € £ 2 (R d ) and 6 G K <i><d an invertible matrix. Then the function 
u g defined in l|22fl satisfies that 

(u g ) b (u) = | dot fo| 1/2 ( X] |%(6(w + s))| 2 ) 1/2 = l a. e.ueE(s,b). 

In particular, if fi(E(g,b) n [0, = 1 i/ien £/ie function ip defined by tp — u g has 
orthonormal (b*)~ 1 Z d translates (i.e. {T^.yi^ip : k € Z d } is orthonormal) . 

The next Lemma states that for each invertible matrix 6 the distance in £ 2 (R d ) of 
two arbitrary functions is bigger than the distance of its Grammians in £ 2 ([0, l] d ). 

Lemma 5.3. Let g,he £ 2 {R d ), b e GL d (R) and u g as in (1221) . Then we have 
11.9 ~ h\\c(s. d ) ^ 11.96 - h b\\c 2 ([o,i] d ) and 

11.9 - u 9ll£ 2 (R d ) = 1.95 - XB(g.b)n[0, l] d ll£ 2 ([0,l] d )- 

Proof. By J23J 

11.9 - h \\%(m<i) = 11(5 - /7 -)fcll£ 2 ([o,i] d ) 



/ \detb\y2\g(b(uj + s)) ~ h(b(iu + s))\ 2 dcj 
J[o,i]« 



> llfffc - M£ 2 ([0,l] d ), 

where in the last inequality we used for the £2^) norm the inequality \\x — y\\ > 

IN -Nil- 

For the second equation, we compute directly 

2 

du> 



JbE(g,b) 



(|det6|Ea l.9(^ + MI 2 ) 1/2 
3HI 2 



/6B(9.b) |deto|X; s \g{uj + bs)\ 2 
If we now pcriodizc and change variables, we obtain the result. □ 



\ 1/2 

det6|^|.g(cj + os)| 2 I -1 



dcu. 
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Remark. Note that if n(E(g,b) f~l [0, l] d ) = 1 then 

11.9 ~ u g\\c 2 {R d ) = ll.9t. — X[o,i] d II £ 2 ([o.i] d ) • 

A consequence of Lcmma l5.3l is the following: Assume that we want to be able to 
find a function ip that is close to a given function f £ Sd and that has orthonormal 
translates with respect to a lattice (b*)~ 1 Z d . Then we will need that /{, is close (in 
£ 2 ([0, l] d )) to the constant function 1. Now, if for a given matrix b, ft is far from 
X[o,i] d i will the choice of a different matrix b improve the error? The next theorem 
establishes the interesting result that fb in fact converges almost everywhere and 
also in norm to X[o,i] d ll/lh when ||6|| — > 0. 

There is a natural interpretation in time domain of this result: when ||6|| becomes 
small, then grows, that is the associated lattice (b*)~ 1 Z d becomes sparser 

and since our functions are in C 2 (M. d ) they will have some decay at infinity which 
will imply that the scalar product between two of its translates will be small. 

In fact the theorem is more general. The convergence also holds for the p- 
Grammians of a function / G C p that we denote by fb,p and are defined as 

f b>p (u)) = (|det6|£ sezd \f{b{uj + s))\P) 1/p . Note that for p = 2, h.z coincides 
with our previous /&. 

The following theorem is a generalization of a result proved in |.TWW05] for the 
C 1 case in a completely different context. 

Theorem 5.1. For any f G C p (R d ), 

fb, P —> ll/llpX[o,i]< a.e. and inC p (\Q,\] d ) when \\b\\ -> 0. (23) 

We postpone the proof of the theorem until the last section. 

Let us now see, how we deduce from Lemma 15.31 immediately that when the 
matrix b is fixed, that is the lattice T = (6*) _1 is fixed, then the set Wo(r) is not 
dense in the sphere: By Parseval, Lemma 15.31 and using that g b = 1 a. c. for every 
g G Wo(r) we have for f e S d and g G W (r): 

11/ - 9\\c 2 (R d ) > Wfb - 9b\\c 2 ([0,l] d ) = Wfb - X[0,l] d ll£ 2 ([0,l] d ) 

The proof is completed choosing a function / in Sd whose Grammian (with respect 
to b) is far from X[o.i] d - 

The following argument from Yang Wang [private communication] shows that 
if the dilation a is fixed, then Wo(o-) is not dense in Sd- Assume that Wo (a) is 
dense in Sd. Let / G Sd be an arbitrary function such that < /, D a f >^ 0. Choose 
ipn £ Wo (a) such that ip n — > /, n — > 00 in Sd ■ Then we have: 

= < Vn, D a 1p n >-►</, D a f >, 

which is a contradiction. 

Now we are ready to prove a density result for the set of generators of orthonor- 
mal (not necessarily complete) wavelet systems: 

Theorem 5.2. The set Wo of generators of orthonormal wavelet systems is dense 
in Sd in the induced £2(R d ) metric. 

Proof. Let / G Sd and e > be given. Choose r, R, g and h as in Theorem 13.21 
By Theorem 15. II there exists b G GLd(M.) with small enough norm, such that 

Wb - X[o,i] d \\c 2 {[o,i] d ) < 9 
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and hb > a.e. Note that this is possible since \h\ > in {uj : r < \\u\\ < R}. Now 
using Lemma 15.31 

s 

\\h - u h \\ C 2 {Rd) < - 

and therefore the inverse Fourier transform of Uh has orthonormal translates in the 
lattice (b*y 1 Z d . 

By choosing the dilation a to be a = — 7<j, the set 

{\dctay/ 2 u h {a 3 uj)e- 2lTa3u)(b ' rlk :keZ d ,j£ z} 
is orthonormal, and consequently, if ip is such that ip = Uh 

U det a\ j / 2 ip(a j x - (6*) _1 fc) : k £ Z d , j e z| 
is an orthonormal wavelet system, and so ip £ Wo- Further we have 

II/- ^|| < \\f-g\\ + \\g-h\\ + \\h-u h \\ <e 
which shows that Wo is dense in £ 2 (R d ). □ 
5.1. Proof of Theorem 15. U 

Proof. We will divide the proof in several steps. Let us denote the unit cube by 
Q = [0, l] d . Wc will denote by || • || p the p-norm in C p (M. d ), and || ■ \\cp(q) the p-norm 
in &(Q). 

• We first prove the theorem for the case that / = \i where 7 is a finite d- 
dimcnsional interval. In this case, we compute 

f b {u) = I dct b\ Xi{b(u + s)) = \ dct b\ + «)• (24) 

S S 

Let 

N^b) := {s £ Z d :bQ + bs C /} 

N a {b) := {s 6 Z d : fj,((bQ + bs) n 7) > and ^((6Q + 6s) n I c ) > 0}, 
and call rij = and n Q = #^0. Then, 

X! X6-i/(w + s) = (rij + m(6, u>)), m(b, w)eN,0< m(6, w) < n Q . (25) 

Note that | dot 6| m(b) < < | det6|(n,-(&) + n (b))\ and hence 

< fi(I) - I det b\ rn(b) < I det b\ n a (b). (26) 

We will now see that | det 6| n (i>) — » for ||6|| — > 0. For this take e > 0, and 
let 6 6 GX d (M) be such that ||6|| < er/diam(Q). Since ||6(a; - y)\\ < \\b\\\\x - y\\ 
then diam(feQ) < ||&||diam(Q) < e. Therefore, if s £ N (b), bQ + bs C {uj £ M d : 
dist(w,<97) < e}. Now 

I det &| no (6) = + 6s ) - e Rd : dist ( w > dI ) < £ D < ec i- ( 27 ) 

sGJV„ 

Therefore we see from l|26[) that 

|det&|n;(&) ->/i(7) when ||6|| -> 0, (28) 
and using (|26() wc have \dctb\J2 s Xb- 1 i{ Ll -> + s ) —* M-O a.e. Furthermore, fbp 

Ki) 1/P = \\f\\ P - 
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In addition, since det b\ J2 S Xfc- 1 /^ + s ) < | det 6|(ni(6) + n {b)) < C for all b 
such that ll&H < 5 < 1 we can apply dominated convergence to obtain 



f btP - Ki) 1/p x Q 



when ||b|| -> 0. 



(29) 



• Now we will prove the theorem for the case that / is a finite linear combination 
of characteristic functions of intervals, i.e. / = X)j=i a jX/j i where Ij are almost 
disjoint (i-dimensional intervals (here by almost, we mean that the intersection can 
have at most measure 0). By the disjointness of the Ij we immediately obtain 

and by the previous item, this converges to ||/|| p . 

In addition, because of the uniform boundcdncss of fb iP with respect to &, we 
obtain convergence in C P (Q) by the dominated convergence theorem. 

• Assume now that / > 0, / <E C p (R d ). Let e > and < g n /" / a.e., with g n as 
in step 2. 

|/6,» - ll/IU < 

|/&,pM - (g n )b,p( W )\ + KflO&.pM - ll.9n|lpl + lllffnllp ~ ll/llpl (31) 

First note that, as in Lemma I5"3l for any f,g S £ p (M d ) we have ||/6, p — gb,p\\c(Q) ^ 
||/ — ,g and therefore 

/ \ IMUp(q) , /qo ^ 

K9n)b,p '' Jb,p- \<J Z ) 

We can therefore choose a subsequence g„. , such that (g nj )b, P (to) — > fb, P (u>) a.e. in 
<5 when j — > oo. Hence, for a large enough j, the first and third term in (|31fl are 
each less than e/4. 

Now, for a fixed j, we can choose <5 such that, if ||6|| < 5, then \{g nj )b,p{<^) — 
WdnjWpl < £/2 which gives the point-wise convergence. 

For the convergence in £ p (M d ), using the previous results, we write 

11/6* - ll/UU(Q) < 

\\fb,p ~ (9n)b,p\\cP(Q) + \\(9n)b,p ~ \\9n\\p\\ CP{Q) + \\\\9n\\p - ll/llpll £P (Q) ■ (33) 

IMIp 

By the particular choice of the sequence g n , we have that g n > /. Therefore the 

first and last term of l|33|) go to when n — > oo. The middle term goes to by the 
previous step. 

Since fb. p = \f\b,p the result holds for arbitrary / S C p (M. d ). □ 

5.2. Hilbert Spaces. The results in this paper carry over to abstract separable 
Hilbcrt spaces, via unitary isomorphisms, where the translation and dilation oper- 
ators are replaced by arbitrary unitary operators. General systems obtained in this 
context (unitary systems) have been studied in detail in |DL98) . 

Finally, we would like to mention that M. Bownik |Bow05| . independently of 
our work and using different techniques has obtained for the case of wavelet frames 
some similar density results as in section 3. 
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